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Abstract 

We perform a comprehensive study of the Higgs sector of the Minimal Supersymmetric 
extension of the Standard Model in the case where all Higgs bosons are rather light, 
with masses of 0(100 GeV), and couple maximally to electroweak gauge bosons and 
strongly to standard third generation fermions, i.e. for large values of the ratio of the 
vacuum expectation values of the two Higgs doublet fields, tan f3. We first summarize 
the main phenomenological features of this "intense-coupling" scenario and discuss the 
available constraints from direct searches of Higgs bosons at LEP2 and the Tevatron as 
well as the indirect constraints from precision measurements such as the p parameter, 
the Zbb vertex, the muon (g — 2) and the decay b — ► sy. We then analyze the decay 
branching fractions of the neutral Higgs particles in this regime and their production 
cross sections at the upcoming colliders, the Tevatron Run II, the LHC, a 500 GeV 
e + e~ linear collider (in the e + e~ and 77 options) as well as at a fA + /j,~ collider. 
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1. Introduction 



A firm prediction of the Minimal Supersymmetric Standard Model (MSSM) [1] is that, among 
the five scalar particles which are present in the extended Higgs sector [2], i.e. two CP-even 
h and H, a pseudoscalar A and two charged if ± bosons, at least the lightest Higgs boson h 
must have a mass below 130 GeV when radiative corrections are taken into account [3-5]. 
If the minimal version of Supersymmetry is indeed realized in Nature, this particle should 
be therefore accessible at the next generation of high-energy colliders, the high-luminosity 
Tevatron [6], the LHC [7, 8] and a future e + e~ linear collider [9, 10]. 

In the decoupling regime [11], that is when the Higgs bosons H , A and are very 
heavy implying that they are almost degenerate in mass, Ma ~ M H ~ M H ±, the lightest 
Higgs particle h will have the properties of the Standard Model (SM) Higgs boson. The 
experimental search for this particle will therefore be straightforward and there is now little 
doubt, in view of the detailed phenomenological and experimental analyses performed in 
the recent years, that it will not escape detection and that at least some of its fundamental 
properties can be pinned down. Unfortunately, the other MSSM Higgs bosons will have 
masses of 0(1 TeV), and will therefore be too heavy to be accessible directly, implying that 
one would not be able to distinguish between the MSSM and the SM Higgs sectors. 

An opposite and much more interesting situation would be the one where the mass of 
the pseudoscalar A boson is not much larger than the maximal value allowed for M^. In 
this case, the three neutral Higgs bosons h, H and A and the charged Higgs particles will 
have comparable masses, Mh ~ Mh ~ Ma ^ 130 GeV and M^± < 150 GeV. A particularly 
interesting scenario is when the parameter tan /3, the ratio of the vacuum expectation values 
of the two Higgs doublet fields which break the electroweak symmetry in the MSSM, is 
large 1 , leading to a CP-even Higgs boson $^ [either h or H] which is almost degenerate 
in mass with the pseudoscalar A boson. As a consequence of the SUSY constraints on the 
Higgs sector, this CP-even $^4 boson will have almost the same couplings as the A boson, 
and therefore will couple strongly to the third generation 6-quarks and r leptons [since the 
couplings are proportional to tan/3], while the couplings of the Z-bosons to $aA pairs [and 
the corresponding couplings of the W bosons to <&aH ± pairs] will be maximal, i.e. not 
suppressed by (sine or cosine of) mixing angle factors. The other CP-even Higgs particle, 
that we will denote by will have almost the couplings of the SM Higgs boson, i.e. its 
couplings to weak gauge bosons and to top quarks are not strongly suppressed by mixing 
angle factors and are therefore almost maximal. 

1 Values tan (3 > 3-10, depending on the mixing in the scalar top sector, are required to maximize the 
h boson mass and to evade the experimental constraint from LEP2 searches [12-14], as will be discussed 
later. Very large values of tan/3 <~ m t /mf, <~ O(50) are favored if one requires Yukawa coupling unification 
at the GUT scale; see Ref. [15]. In the constrained MSSM or minimal Supergravity model [16], large values 
of tan/3 lead naturally to rather light A, H and bosons; see e.g. Ref. [17]. 
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This scenario, with all MSSM Higgs bosons being light with almost maximal couplings 
to gauge bosons and strong couplings to third generation fermions, will be called hereafter, 
the "intense-coupling" regime. 

In this non-decoupling scenario, all three neutral Higgs bosons, as well as the charged 
Higgs particle, would be accessible at the next generation of experiments. The experimental 
searches for the neutral Higgs bosons in this case will be slightly more involved and much 
more interesting than in the decoupling regime, since a plethora of production and decay 
processes, with rates which can be violently different from the ones of the SM Higgs particles, 
have to be considered and studied in detail in order to detect individually all the three 
particles and to determine their basic properties. In addition, two features might render 
the situation somewhat more delicate: the h, H, A boson masses can be very close to each 
other, implying that backgrounds for a given Higgs boson signal would come from another 
Higgs signal and the total decay widths of some of the Higgs bosons can be rather large [in 
particular for very large values of tan j3\ implying broader signals. 

The search for the charged Higgs particle would be straightforward. Indeed, since in the 
MSSM the A and ^ boson masses are related by M 2 H± ~ M\ + M^, M H ± will be smaller 
than ~ 150 GeV for Ma ^ 130 GeV, implying that if ± can always be produced in top quark 
decays, t — > H + b, and can easily be detected at hadron [6, 7] or e + e~ colliders [9, 10]. 

The purpose of this paper is to perform a comprehensive study of the MSSM Higgs 
sector in the intense-coupling regime and the implications of the near mass degeneracy 
M h ~ M H ~ M A , for the search of the Higgs bosons at the Tevatron, the LHC, a 500 GeV 
e+e~ collider and a collider 2 . We will first discuss the parameterization of the Higgs 

sector and derive rather simple, accurate and useful expressions for the radiative corrections 
to the Higgs masses and couplings which are valid in the region 90 GeV < M A < 130 
GeV for moderate (~ 10) and large (~ 30-50) tan (3 values. We will then analyze in this 
"intense-coupling" regime, the various decay branching ratios of the h, H and A particles 
and collect all their production cross sections at these machines. We will not only analyze 
the main production channels which allow to detect these particles, but also sub-leading 
processes which would allow for the measurement of some of their fundamental properties. 
The various experimental signatures for the Higgs particles will be discussed. 

Before that, we will discuss in detail the various constraints on this intense-coupling 
regime with high tan (3 values: the direct constraints from the Tevatron and LEP2 searches 
for MSSM Higgs bosons [and the implications of the possible ~ la evidence [12] for a SM- 
like Higgs boson with a mass around 115 GeV] as well as all the indirect constraints from 

2 Several articles in the literature have dealt with some aspects of the scenario discussed here, and in 
particular, with the consequences of large tan/3 values for MSSM Higgs boson production at future colliders 
[18, 19]. The conclusions of these papers overlap with the ones we obtain here. 
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precision measurements [20, 21], i.e. the W boson mass, the effective electroweak mixing 
angle sin 2 6 W and the Z boson decays into bb final states, the radiative decay b — > S7 [22] 
and the anomalous magnetic moment of the muon, (g — 2)^ [23]. 

The paper is organized as follows. In the next section, we will present the physical set- 
up and summarize our parameterization of the MSSM Higgs sector, and display the Higgs 
boson masses and couplings in the intense-coupling regime. In section 3, we discuss the 
implications of the LEP2 and Tevatron direct Higgs searches and the indirect constraints 
from high precision measurements on this scenario. The various decay branching fractions 
of the neutral Higgs bosons will be analyzed in section 4. The production cross sections of 
the neutral Higgs bosons at the Tevatron, the LHC, a 500 GeV e + e~ linear collider and a 
fi + fi~ collider will be given in section 5, together with a discussion of the various processes 
and signatures which can be used. A conclusion will be given in section 6. 

2. The MSSM Higgs sector and the intense— coupling regime 
2.1 Radiative corrections in the Higgs sector 

In the MSSM, besides the four Higgs boson masses, two mixing angles define the properties 
of the scalar particles and their interactions with gauge bosons and fermions: the ratio of 
the two vacuum expectation values tan/5 = v 2 /vi of the Higgs doublets and a mixing angle 
a in the neutral CP-even sector. Supersymmetry leads to several relations among these 
parameters and, in fact, only two of them, taken for convenience to be Ma and tan/3, are 
independent. These relations impose, at the tree-level, a strong hierarchical structure on the 
mass spectrum [e.g. Mh < Mz] some of which are, however, broken by radiative corrections. 
The leading part of these corrections grows as and logarithmically with the common stop 
mass Ms [3]. Including the mixing in the stop and sbottom sectors and all the subleading 
terms, the full expressions of the radiative corrections become rather involved [4, 5]. 

To illustrate the qualitative behavior of the Higgs boson masses and couplings, we will use 
simple analytic expressions [comparable in simplicity to those obtained in [3] for the leading 
m\ piece] where the stop mixing and some of the subleading terms are incorporated. In 
many cases, these expressions give rather good approximations for the neutral Higgs boson 
masses [at the percent level] and couplings [at the ten percent level] compared to what is 
obtained with the full set of corrections in the renormalization group improved approach of 
Ref. [4]. This is sufficient for our purpose in the present section, since our aim is simply 
to understand qualitatively the phenomenological implications of almost degenerate Higgs 
bosons with enhanced Yukawa couplings. 

For our numerical analysis throughout the paper, we will use the program HDECAY [24], 
which calculates the Higgs spectrum and decay widths in the MSSM. The complete radia- 
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tive corrections due to top/bottom quark and squark loops within the effective potential 
approach, leading NLO QCD corrections [through renormalization group improvement] and 
the full mixing in the stop and sbottom sectors are incorporated using the analytical ex- 
pressions of Ref. [4]. [Note also that we will include the leading one-loop and two-loop 
corrections to the Higgs self-couplings.] However, even in this case, our treatment of the 
radiative corrections in the Higgs sector will not be complete, since two sets of potentially 
large corrections are presently not included in the program: 

- The SUSY-QCD and dominant Yukawa radiative corrections to the bottom and top 
quark Yukawa couplings which, in many cases, can be rather important. In particular, 
the corrections to the bottom Yukawa couplings can be very large for high values of tan j3 
and the higgsino mass parameter fi [25]. The net effect of these corrections can be viewed 
as an increase or decrease of the bottom quark mass, depending on the sign of /i. The 
phenomenological impact of these corrections has been already discussed in the literature, 
see e.g. Ref. [18, 25], and we have little to add. 

- Recently, the two-loop 0(af) and 0{a s a^) Yukawa corrections to the Higgs boson 
masses have been calculated [26]. The former set of corrections can increase the lightest h 
boson mass by several GeV in some areas of the MSSM parameter space, while the later 
corrections can be sizeable for large tan (3 values. 

We have nevertheless verified that the numerical results that we obtain for the MSSM 
Higgs spectrum, and in particular for the lightest h boson mass, are rather close to those 
obtained from the complete results of the Feynman diagrammatic approach implemented 
in the program FeynHiggs [27], which include some of the Yukawa corrections of Ref. [26]. 
The difference, of the order of a few GeV, is of the same size (within a factor of two) as 
the expected theoretical error on the determination of the Higgs boson masses. However, as 
already pointed out, this difference will not be of utmost importance for our analysis, and 
our general conclusions will not be altered in a significant way. 



2.2 The Higgs boson masses 



Before discussing the Higgs boson masses in the intense-coupling regime, i.e. when the A 
boson is rather light and tan (5 is very large, let us first present the rather simple analytical 
expression which approximate the Higgs boson masses. We first define the quantity e, which 
parameterizes the main radiative correction, 
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where m t is the running MS top quark mass to account for the leading QCD corrections, A t 
is the stop trilinear coupling, and as the strong coupling constant. The Higgs boson masses 
can be approximated to an accuracy of the order of a few percent compared to the complete 
result [in particular in the case where the splitting between the two stop masses, and to a 
lesser extent sbottom masses, is not very large], as functions of M A , tan/5 and e with the 
expression 
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This radiative correction pushes the maximum value of the lightest h boson mass upwards 
from Mz by several tens of GeV: in the so-called maximal mixing scenario A t = A t — 
fj,/ tan/5 = v^SMs and with M A and Ms of about 1 TeV, one obtains an upper mass bound, 
Mh < 130 GeV for tan/5 ^> 1, a bound that is comparable to the one obtained including the 
full set of corrections in the RG improved approach [4]. 

In the case of the mixing angle a of the CP-even Higgs sector, the correction eq. (1) is 
not accurate enough and one needs to introduce another correction which involves the ratio 
n/Ms, where \x is the higgsino mass parameter. Defining the parameter e' to be 
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the mixing angle a is given in terms of Ma, tan (3, e and e', by 
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Here again, the accuracy of the formula, compared to the case where the full corrections are 
taken into account, is at the percent level if the splitting between the stop (and to a lesser 
extent, sbottom) masses is not too large. 

[In the case of the charged Higgs boson mass, one can also derive a very simple expression 
for the radiative corrections which gives a result that is accurate at the percent level, 



M H ± = \JM\ + M^ - e + with e + = 



3G F M^t 
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where m& is the running MS b mass at the scale of the top quark mass.] 

Let us now discuss in some details the neutral Higgs boson masses in the intense-coupling 
regime, taking into account the dominant radiative corrections e and e'. In fact, this regime 
can be defined as the one where the two CP-even Higgs bosons h and H are almost degenerate 
in mass, M h ~ M H . We will therefore first concentrate on the Higgs sector in the limit 
M h = M H . 
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Solving the equation (2) for M\ — M% = 0, which is a second order polynomial equation 
in the variable M\, 

M\ + 2M 2 A [M 2 Z {1 - 2 cos 2 2(3) + e cos 2(3} + M| + e 2 - 2Mfecos2/3 = (6) 

one obtains a discriminant A' = — sin 2 2(3{2M Z cos 2/3 — e) 2 < 0. The only way for the 
solution to be real is therefore to have either sin 2(3 = or e = 2Mf cos2/3. The last 
possibility gives M\ = —M\ which has to be rejected, while the former possibility gives 
M\ = Mf + e with (3 — | . In fact, this solution or critical mass corresponds to the maximal 
value allowed for Mh and the minimal value that Mh can take, 

M c = AC ax = Mf n = \j M z + e (7) 

In addition, in the large tan (3 regime, eq. (2) for the masses of the CP-even Higgs bosons 
simplifies to 

Kh = \{Ml + Ml + eT\M 2 A -M 2 z -e\) (8) 

which means that 

M A > M c =>• M H = M A and M h = M c 

M A < M c =>- M h = M A and M H = M c (9) 

and therefore the A boson is always degenerate in mass with one of the CP-even Higgs 
bosons, and if the masses of the latter are equal, one has M H = M h = M A = M c . 

2.3 The Higgs boson couplings 

We turn now to the couplings of the Higgs bosons, which determine to a large extent, the 
production cross sections and the decay widths. The pseudoscalar Higgs boson has couplings 
to isospin down (up) type fermions that are (inversely) proportional to tan (3 and, because of 
CP-invariance, it has no tree-level couplings to gauge bosons. In the case of the CP-even h 
and H bosons, the mixing angle a enters in addition. The couplings of the CP-even Higgs 
bosons to down (up) type fermions are enhanced (suppressed) compared to the SM Higgs 
couplings for values tan/9 > 1; the couplings to gauge bosons are suppressed by sin(/5 — a) 
or cos(/5 — a) factors; see Table 1. In fact, the h and H bosons share the couplings of the 
SM Higgs to the gauge bosons; also, the squared sum of the couplings to fermions do not 
depend on a and obey the sum rules: 

9hM + 9Hdd = V cos2 P » 9luu + 9 2 huu = V sin 2 (3 , g 2 hVV + g 2 HVV = 1 (10) 
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Developing the formulae for sin 2a and cos 2a for large values of tan/5, tan (3 3> 1, 

• « • ^ M l + M i- e 7 sin2 /3 n nn M\ -Ml + e/cos2/5 
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The couplings of the neutral Higgs bosons are displayed in Table 1 together with their 
limiting values for tan/5 ^> 1 and M A > Mc (upper values) or Ma < Mc (lower values). 
The accuracy is of the order of 10% for A t ,At> and [i of the order of 1 TeV and improves 
in the case of zero mixing. One can see that, not only in the decoupling limit Ma — ► 00 
but also in the intense-coupling regime with tan/5 ^> 1, the CP-even Higgs boson with a 
mass Mc has, up to a sign, SM-like Higgs couplings to gauge bosons and up-type quarks. 
For down-type fermions, the situation is more complicated: in general the coupling is larger 
than the SM Higgs coupling, but in some cases, when all terms add up to zero [i.e. when 
fi ~ /itan/5 in our approximation], the coupling becomes very small and even vanishes. 
The other Higgs boson, with a mass M A) has couplings of the same order as the A boson, 
i.e. it has rather small couplings to gauge bosons and couples to up-type and down-type 
fermions proportionally to, respectively, =Fl/ tan/3 and tan/5 [in the case of small 

In the intense-coupling regime, we will therefore define a SM-like Higgs boson § H and 
a pseudoscalar-like Higgs boson $^ as follows: 



^ H = h if M A > M c and <$> H = H if Ma < M c 
^ A = h if M A < M c and $ A = H if Ma > M c 



(13) 





9&UU 


9*dd 




h 


cos a/ sin/5 Z- fl+ f 2 /t^p 


-sin a/ cos /5 ^J n/3 


sin(/5-a) Z-f 1+{ i +f2) /t,n3 


H 


1 • — *fi — H 1 tan 3 

sin a:/ sm/5 ' 


1 n — >tan 8 

cos a/ cos (3 _^_ h ^ np+h 


cos(/5-a) Zt\ +{1 - h)/tanP 


A 


l/tan/5 


tan/5 






Table 1: Neutral Higgs boson couplings to fermions and gauge bosons in the MSSM normal- 
ized to the SM Higgs boson couplings, and their limits for tan (3 >> 1 for Ma > M c ( upper 
values) and Ma < M c (lower values). 
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An additional set of couplings that we will need in this study are the trilinear couplings 
among the neutral Higgs bosons. These couplings will be defined in units of A = M\jv and 
are renormalized not only indirectly by the renormalization of the angle a but also directly 
by additive terms proportional to the parameter e. In the e approximation [which here, gives 
only the magnitude of the Higgs self-couplings, i.e. a few ten percent in general but in some 
cases, more than a factor of two] and keeping only 1/ tan/5 terms, one has: 
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(14) 



Again, the upper (lower) values are the limits for tan/3 ^> 1 with larger (smaller) 
than Mc- The SM trilinear self-coupling 3Mj^ /M§ is approached in the case of h in the 
decoupling limit, and in the case of H (up to a sign) in the intense-coupling limit. Without 
radiative corrections, one of these Higgs bosons would have A-like couplings while the other 
would decouple for large tan /3 values. 

If 90 GeV < M4 < 130 GeV, one would have also a rather light H ± boson with a mass 
120 GeV < M H ± < 150 GeV. Its coupling to fermions is a parity violating mixture of scalar 
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and pseudoscalar currents, gn+ud oc (1 — 7s)m M / tan/3 + (1 + 75)m d tan/9. As in the case of 
the A boson, the couplings to down (up) type fermions are (inversely) proportional to tan (3. 
Only the couplings to the (i, b) and (r, v T ) isodoublets are relevant in most cases. We will 
not discuss the properties of this particle further in this paper. 

Finally, note that the couplings of the CP-even h and H bosons to ZA and W + H~ pairs 
are proportional to cos(/3 — a) and sin(/3 — a), respectively, while the W + H~ A coupling is 
not suppressed by these factors. This means again, that for large tan [3 values the CP-even 
$a boson with a mass close to M A couples maximaly to ZA and W Zf H ± states while the 
SM-like boson has small couplings. 

2.4 Inputs for the numerical analysis 

As mentioned previously, our numerical analyses all through this paper will be based on the 
program HDECAY for the evaluation of the masses, couplings and decay branching ratios. We 
will use the running mass of the pseudoscalar Higgs boson A, together with tan f3 and the 
relevant MSSM parameters which enter in the radiative corrections [A t , \x and the soft SUSY 
breaking third generation squark masses, = iTijj = iTi£> = Ms] as inputs, to obtain 
the masses of the two CP-even h and H bosons and the charged particles. The pole 
bottom and top quark masses will be fixed to m b = 4.25 GeV and m t = 175 GeV. In most 
of the analysis, we will use two values of tan (3, a moderate and a large one, tan (3 — 10 and 
30, and deal with the maximal mixing scenario with A t ~ y/6M$ with Ms fixed to 1 TeV. 
For the Higgs boson masses however, we will also discuss the no-mixing and typical mixing 
scenarii with, respectively, A t = and A t = Ms, and give some illustrations for a small and 
very large value of tan (3, 5 and 50. 

Because of the approximations discussed in section 2.1 [in particular because the SUSY 
radiative corrections to the Yukawa couplings are not incorporated in the present version of 
HDECAY] the other SUSY parameters play only a minor role and we will fix them to A b — 1 
TeV and /j, = M 2 = 350 GeV for the discussion of the Higgs boson masses, couplings and 
decay branching ratios and A b = /x = M 2 = 1 TeV for the production cross sections. We have 
verified that within our approximation, this difference does not lead to significant changes 3 . 

The neutral Higgs boson masses that we obtain are shown in Fig. 1 for a mass of the 
A boson varying from 90 to 130 GeV for four representative values of tan/3 from small to 
very large, tan/3 = 5,10,30 and 50, in the three scenarii of stop mixing discussed above. 
The normalized Higgs boson couplings to fermions and gauge bosons are shown in Fig. 2 for 
tan (3 = 10, 30 and 50 in the scenario with large mixing, A t = 2.6 TeV. 

3 However, we recall again, that if the SUSY corrections to the &-quark Yukawa coupling are included, this 
might not be the case for large tan/3 values [18, 25]. This is particularly true in the discussion of the branching 
ratios: while in our case, BR($ — > 66)/BR($ — > t+t~) for instance is simply given by ~ 3m^(M|)| SM /m^, 
it can increase or decrease for large values of tan/3 depending on the sign of zz. 
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Figure 1: The masses of the lightest and heavier CP-even Higgs bosons h and H [in GeV] 
as a function of Ma for four values oftanf3 = 5, 10,30 and 50 in the case of maximal mixing 
(upper), typical mixing (central) and no-mixing (lower). 
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Figure 2: The couplings of the CP-even Higgs bosons h and H [normalised to the SM Higgs 
couplings] as a function of the masses for three values of tanf3 = 10, 30 and 50, in the 
maximal mixing scenario. 
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3. Experimental Constraints on the intense— coupling scenario 
3.1 Direct constraints from LEP2 and Tevatron searches 

The search for the Higgs bosons was the main motivation for extending the LEP2 energy 
up to y/s ~ 209 GeV [12]. In the SM, a lower bound M H o > 114.1 GeV has been set at 
the 95% confidence level, by investigating the Higgs-strahlung process, e + e _ — > ZH° [13]. 
In the MSSM, this bound is valid for the lightest CP-even Higgs particle h if its coupling 
to the Z boson is SM-like, i.e. if 9zzh/ 9zzh° = s ^ n2 (P — a) — 1 [almost in the decoupling 
regime] or in the less likely case of the heavier H particle if 9zzh/9zzh° = cos2 (P — a) ~ 1 
[i.e. in the non-decoupling regime with a rather light M H \\ see e.g. Ref. [28]. 

A complementary information is obtained from the search of Higgs bosons in the asso- 
ciated production processes e + e~ — > Ah where the 95% confidence level limits, M h > 91.0 
GeV and Ma > 91.9 GeV [14], on the h and A masses have been set. This bound is ob- 
tained in the limit where the coupling of the Z boson to hA pairs is maximal, g 2 ZhA / 'g Z zH° = 
cos 2 (/3 — a) ~ 1, i.e. in the non-decoupling regime for large values of tan/3. This limit is 
lower than the one from the Higgs-strahlung process, due to the less distinctive signal and 
the P 3 suppression near threshold for spin-zero particle pair production. Note that for small 
Ma and large tan/3 values, M H becomes small enough, cf. Fig. 1, in the no-mixing scenario 
to allow for the possibility of the process e + e~ — > HA, which is suppressed by sin 2 (/3 — a). 

Deriving a precise bound on M h for arbitrary values of Ma and tan (3 [i.e. not only in the 
decoupling limit or for tan (5 >> 1] and hence, for all values of the angle a, is more complicated 
since one has to combine results from two different production channels, which have different 
kinematical behavior, cross sections, backgrounds, etc.. However, some exclusion plots for 
sin 2 (/3 — a) versus M^ from the Higgs-strahlung process [and which can be used to constrain 
the mass of the H boson if sin 2 (/3 — a) is replaced by cos 2 (/3 — a)] and cos 2 (/3 — a) versus 
Ma + Mh [with Mh ~ M A ] from the pair production process, have been given in Refs. [13] 
and [14], respectively 4 . We have fitted these exclusion contours and delineated the regions 
allowed by the LEP2 data up to ^fs = 209 GeV in the [M^tan/?] and [M^, tan/3] planes in 
the no-mixing, typical mixing and maximal mixing scenarii [the other SUSY parameters are 
fixed as in section 2]. The domains allowed by the LEP2 constraints 5 are shown in Fig. 3. 

4 The exclusion plot cos 2 (/3 — a) versus Mh + Ma has been obtained with the assumption that the total 
widths of the Higgs bosons are rather small, which is not the case for large values of tan/3 > 30. We will 
assume here, that one can nevertheless extrapolate these limits up to values tan/3 ~ 50. 

5 These domains agree rather well with those derived by the LEP collaborations in their combined analyses 
[13, 14]. This means that although we use different programs for the implementation of the radiative 
corrections [HDECAY with an improved RG Higgs potential versus a complete implementation of the two- 
loop radiative corrections in the Feynamn diagramatic approach] , the obtained Higgs boson spectrum is the 
same with a very good approximation, 0(2 GeV). This also means that our fitting procedure of the LEP2 
contraints on the Higgs boson masses versus the mixing angle is rather good. 
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Figure 3: The allowed regions for Ma [left] and M h [right] from LEP2 searches as a function 
of tan (3 (colored regions, or black, dark grey and light grey regions) in the case of maximal 
(bottom), typical (central) and no-mixing (up). The red (or dark grey) regions indicate where 
114 GeV< Mh < 117 GeV and sin 2 (/3— a) > 0.9 and the green (or light grey) regions indicate 
where 114 GeV< M H < 117 GeV and cos 2 (/3 - a) > 0.9. 
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In these figures the colored (blue, green and red) areas correspond to the domains of the 
parameter space which are still allowed by the LEP2 searches. As can be seen, the allowed 
regions are different for the three scenarii of mixing in the stop sector. 

In the no-mixing case, values tan (3 < 10 are excluded for any value of the pseudoscalar 
mass M A > 90 GeV. For a light A boson, M A < 100 GeV and tan (3 > 10, the masses 
M h ~ M A pass the bound ( > 91 GeV) from the Higgs pair production processs, while H is 
light enough (Mh < 115 GeV) to be probed in the Higgs-strahlung process. For increasing 
M A , this situation holds only for larger tan/3: for instance, for M A > 115 GeV, one needs 
values tan (3 > 40 to cope with the experimental bounds [i.e. to make H light enough to be 
produced]. In the typical mixing case, values tan/5 < 10 are excluded for M A close to 90 or 
130 GeV. For intermediate M A values, one needs larger tan f3 values to evade the constraint 
from the Higgs-strahlung process since the ghzz coupling is increasing while Mh is relatively 
small, leading to a sizeable a(e + e~ — > hZ). [Note that here, there is a strong interplay 
between the variation of sin 2 (/5 — a) and Mh in the phase space in cr(e + e~ — > hZ), which 
explains the bump around M A ~ 115 GeV.] For large M A , tan/5 values close to 5 are still 
allowed since M h can be larger than 114 GeV. Finally, in the maximal mixing case, values 
tan (3 > 7 are excluded for 90 GeV < M A < 130 GeV. However, for large M A , only values 
tan (3 > 3 are excluded since here, the large stop mixing increases the maximal value of M h 
to the level where it exceeds the discovery limit, M h > 114 GeV. 

We turn now to the implications of the 2.1a evidence for a SM-like Higgs boson with a 
mass M^ H = 115.6 GeV, seen by the LEP collaborations in the Higgs-strahlung process [12]. 
In view of the theoretical and experimental uncertainties, this result will be interpreted as 
favoring the mass range 114 GeV < M$ H < 117 GeV. Furthermore, since this Higgs boson 
should be almost SM-like, we will impose the constraint g 2 zz^ H / 9zzh° — 0-9 so that the 
cross section for the Higgs-strahlung process, a(e + e~ — > is maximized. Of course, 

at the same time [and in particular in the case where the "observed" Higgs boson is H] the 
experimental constraints from the pair production process, M A ,M h > 91 GeV, should be 
taken into account. 

The green areas show the regions where the "observed" Higgs boson is the heavier CP- 
even H particle. This occurs only for no mixing with values tan/3 > 15 and M A < 110 
GeV. In this case, Mh ~ 115 GeV and has SM-like couplings, while Mh is still larger than 
91 GeV. For large stop mixing, the H boson is always heavier than 117 GeV and cannot 
be observed. The red areas correspond to the regions where the "observed" Higgs boson is 
in turn the lighter h particle, while Mh is heavier than 117 GeV. These regions are larger 
for the typical mixing scenario: in the no-mixing case, it is unlikely that the h boson mass 
exceeds 114 GeV except for very large tan/5, while in the maximal mixing scenario, the h 
boson is usually heavier than 117 GeV, except for very low tan/5. 
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Finally, let us briefly discuss the constraints on the intense-coupling scenario from the 
Higgs boson searches at the Tevatron. The first important constraint is coming from decays 
of the top quarks into 6-quarks and charged Higgs bosons, t — > bH + . For tan (3 values around 
unity [which are ruled out as discussed above] and for large values of tan /3, the H~tb coupling, 
gn-tb m f / tan/5 + mb tan /3, is very strong. The branching ratio BR(t — > bH + ) becomes 
then comparable and even larger than the branching ratio of the standard mode, t — > bW + , 
which allows the detection of top quarks at the Tevatron. A search for this top decay mode, 
with the H boson decaying into tv t final states, by the CDF and DO Collaborations [29] 
allows to place stringent limits on the value of tan /3 for charged Higgs boson masses below 
~ 150 GeV: for M H ± ~ 120 GeV [i.e. M A ~ 90 GeV] one has tan/? < 50-60, while the bound 
becomes very weak for M#± ~ 150 GeV [i.e M A — 130 GeV] where one has tan/3 < 100 [i.e. 
beyond the limit where the 6-quark Yukawa coupling is perturbative] . 

Furthermore, constraints on the MSSM Higgs sector at the Tevatron can be derived [30] 
by exploiting the data on r + r~+ 2 jets used by the CDF Collaboration [31] to place limits 
on third generation leptoquarks. Indeed, in the associated production of the pseudoscalar A 
boson with bb pairs, pp — > qq/gg — > bbA, with the A boson decaying into r + r~ final states, 
the cross section is proportional to tan 2 /3 and can be very large for tan/3 ^> 1. In addition, 
since in this regime, one of the CP-even Higgs bosons is always degenerate in mass with the 
A boson and has couplings to bb [and t + t~] pairs which are also proportional to tan/3, the 
rate for bbr + T~ final states due to Higgs bosons in this type of process is multiplied by a 
factor 2. Since there is no evidence for non-SM contributions in this final state as analyzed 
by the CDF Collaboration, one can place bounds in the plane spanned by tan/3 and M A . 
The upper bound is tan/3 < 80 for M A > 90 GeV, and is weaker for higher values of M A , 
tan (3 < 100 for M A ~ 130 GeV [30]. Therefore, this constraint is less severe than the one 
due to top decays into charged Higgs bosons [29]. 

There are two other processes for MSSM Higgs boson production which could be relevant 
at the Tevatron as will be discussed later. The gluon-gluon fusion mechanism, gg — > $ with 
$ = h, H or A, can have large cross sections in particular for large tan (3 values, but the 
backgrounds in the main decay channel $ — > bb are too large while for the cleaner 77 decays, 
a x BR($ — > 77) is too small. The associated production of the CP-even Higgs bosons h or 
H with a W or Z boson, has cross sections which are too small, in particular if the gauge 
bosons are required to decay leptonically. These processes will therefore be useful only when 
an integrated luminosity of a few femtobarn will be accumulated, i.e. only at the Run II of 
the machine [see section 5]. 

To summarize, the constraints on the Higgs sector from Tevatron data are not very strong 
and values of tan /3 slightly above tan (3 ~ 50 [and below] are still allowed experimentally. 
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3.2 Indirect constraints from precision data 



Indirect constraints on the parameters of the MSSM Higgs sector, in particular on Ma and 
tan (3, come from the high-precision data: the measurements of the p parameter from M w 
and sin 2 6^, the decays Z — > bb, the muon anomalous magnetic moment (g^ — 2) and the 
radiative decay b — > sj. In discussing these constraints, we will neglect the contributions of 
the SUSY particles which are assumed to be heavy [p — M 2 — 350 GeV and M s = 1 TeV 
with A t = V&Ms]. For an analysis of these SUSY contributions in a constrained MSSM, see 
Ref. [17] for instance. For completeness, we will also give the analytical expressions of the 
Higgs boson loop contributions to the observables in the intense-coupling regime. 

a) The p parameter 

Loop contributions of the MSSM Higgs bosons can alter the values of the electroweak 
observables which have been precisely measured at LEP1, the SLC and the Tevatron [20, 21]. 
The dominant contributions to the weak observables, in particular the W boson mass and 
the effective mixing angle s 2 ^ = sin 2 $w, enter via a deviation from unity of the p parameter 
which, in terms of W and Z boson self-energies at zero momentum transfer, is defined by 
p = (1 — Ap) -1 with Ap = U zz (0)/M z — U ww (0)/Myy Precision measurements constrain 
the contribution of New Physics to be Ap NP < 1.1 • 10" 3 at the la level [20]. 

In the intense-coupling regime, the contributions of the MSSM Higgs sector to the p 
parameter is given by [note that the same expression holds exactly in the decoupling regime]: 



^ Higgs _ ~r — 'M/ 



i z / \ ivi w m w 

with the two functions fi and f 2 given by [32] 



(15) 



In c?v — In x In a; 

+ 



-w 



X Cyy^l X) 



,f2(xi,x 2 )= Xl%2 In— + -(xi + x 2 ) (16) 
X\ — x 2 X\ 2 



with clr = 1 - s 2 w . The first contribution, Ap^f gs = -3G , f M^/(8v / 2tt 2 )/ 1 (M^/M|), is 
simply the one of the SM Higgs boson with a mass M H o = M C) which is close to the Higgs 
mass, (9(100 GeV), favored by the global fits of the electroweak data [21]. The second 
contribution, proportional to f 2 , is the genuine MSSM contribution. The contribution of 
loops involving CP-even and CP-odd Higgs boson, ~ — f 2 (M\/M^ /) M\/M^) in the bracket, 
vanishes since f 2 (x,x) = 0. In fact, only loops which involve particles which have a large 
mass splitting contribute significantly to Ap. This is the reason why the genuine MSSM 
contribution, Ap^!f SM = -G F M 2 ,/(4 v ^7r 2 )/ 2 (M 2 r± /M 2 ,, Mi/M 2 ,), shown in Table 2, is 
always extremely small: the mass difference between the and A bosons is not large 
enough. In turn, the SM-like Higgs contribution has a size which is comparable to the 
experimental error for M H o = M c ~ 127 GeV [we use tan (3 = 50]. 
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M A [GeV] 


90 


110 


130 


A Higgs 
^non-SM 


-0.53 • 10~ 4 


-0.38- 10~ 4 


-0.28 • 10~ 4 


A Higgs 


4.7- 10~ 4 


4.7- 10~ 4 


4.7- 10~ 4 



Table 2: The genuine and the SM-like contributions of the MSSM Higgs sector to the p 
parameter in the intense-coupling regime with tan/3 = 50. 



b) The g— 2 of the muon 

The precise measurement of the anomalous magnetic moment of the muon recently per- 
formed at BNL, a M = - 2 = 11659202(20) • 10~ 10 [23] where we have added the statistical 
and systematical errors, can provide very stringent tests of models of New Physics. In the 
MSSM, the Higgs sector will contribute to a M through loops involving the neutral Higgs 
bosons h, H and A with muons and loops involving the charged Higgs bosons H ± with neu- 
trinos. The contributions are sizeable only for large values of tan/5 for which the 
and H + [iv^ couplings are enhanced; for a recent analysis see Ref . [33] . 

Taking into account only the leading, oc tan 2 (3, contributions [i.e. neglecting the contri- 
bution of the SM-like CP-even Higgs boson <& H ] and working in the intense-coupling regime, 
one obtains for the MSSM Higgs sector contribution to a M [again, the expression holds also 
in the decoupling limit in this approximation] 



Higgs 



G F ml 



tan 2 /? 



m,, 

4^ 



(17) 



This contribution is given in Table 3 for tan/5 = 50 for M A = 90, 110 and 130 GeV, and one 
can see that it is far too small compared to the experimental error. [Here, the approximated 
result which is displayed, is about 10 to 20 % smaller or larger than the exact result]. 



M A [GeV] 


90 


110 


130 


n Higgs 


4.5 • 10~ 12 


2.9 • 10~ 12 


2.0 • 10~ 12 



Table 3: Contribution of the MSSM Higgs sector to for tan/? = 50. 



c) The Zbb vertex 

An observable where the MSSM Higgs sector can have sizeable effects is the Z boson decay 
into bb final states. The neutral Higgs particles h, H, A as well as the charged H ± bosons can 
be exchanged in the Zbb vertex [34] and for large values of tan (3, for which the Higgs boson 
couplings to b quarks are strongly enhanced, they can alter significantly the values of the 
partial decay width T(Z — > bb) [or equivalently the ratio Rb = T(Z — > bb)/T(Z — > hadrons) 
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which is measured with a better accuracy] and the forward-backward asymmetry A b FB [or 
the polarization asymmetry A b LR ]. In terms of the left- and right-handed Zff couplings, 
9l/r = If ~~ e f s wi an d neglecting the 6-quark mass for simplicity, they are given by: 



T(Z -> bb) ~ 



GFMl MY + {g R f] 



V2> 



IT 



3 (gj) 2 - (g e R ) 2 {gj? - (g b R ) 2 
1 1: " 4 (gir + ( 9 %y (g b L y + ( 9 b R y 



A b ~ 



(18) 



In the limit tan/3 ^> 1, the MSSM neutral (N) and charged (C) Higgs boson contributions 
to the left- and right-handed bottom quark couplings are given by 



Sg R /L\N 

tg b R \c 



R/L 
2 r 



Sg 

/ gm b tan (3 
V V2M W 
gm b tan /3 X 2 



-g R C (M H ±,m t ,m t ) 



^9r/l\n + Sg R / L \c 
C 2 (m b ,M A ,M A )±g b L/R C 1 (M A ,m b ,m b ) 

■2s 2 w - l)c 2 (m t , M H ±,M H ±) + gldiMj&M, m t ) 

(19) 



Sg b L \c 



where the functions C\ t2 are given in terms of the Passarino-Veltman two- and three-point 
functions [the latter evaluated at q 2 = Mf with m b ~ 0] by: 



C 2 (m 1 ,m 2 , m 2 ) = C 2 A(m 1 ,m 2 , m 2 ) + -B 1 (m 1 ,m 2 ) 

C 1 (m 1 ,m 2 ,m 2 ) = ~ + 2C 2i (mi, m 2 , m 2 ) - M|C 23 (mi, m 2 , m 2 ) + B^mi, m 2 ) (20) 

The latest experimental values of R b and the forward backward asymmetry are [21]: R b = 
0.21653 ± 0.00069 and A b FB = 0.099 ± 0.002 [note that they deviate by, respectively, +1.1(7 
and —2.4cr from the predicted values in the SM]. This means that the virtual effects of the 
MSSM Higgs bosons should be, in relative size, of the order of 0.3% in R b and 2% in A b FB to 
be detectable. This is far to be the case even for tan (3 values close to 50 and for M A > 90 
GeV as can be seen from Table 4 where the approximate contributions [which are very close, 
at most a few percent, to the exact contributions] are displayed. 



M A [GeV] 


90 


110 


130 


AR b /R b 


-5.0 • 10" 4 


-7.5 • 10~ 4 


-8.3 • 10~ 4 


A^fbMfb 


2.5 • 10~ 3 


2.5- 10" 3 


2.4- 10~ 3 



Table 4: Relative contribution of the MSSM Higgs sector to the two observables R b = T(Z 
bb)/T(Z -> hadrons) and Ap B for tan/5 = 50. 
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d) The b — > 57 decay 

Finally, in the radiative and flavor changing decay b — > 57, in addition to the SM con- 
tribution built-up by W boson and top quark loops, loops of charged Higgs bosons and top 
quarks can significantly contribute in the MSSM together with SUSY particle loops [35, 36]. 
Since SM and MSSM Higgs contributions appear at the same order of perturbation theory, 
the measurement of the inclusive branching ratio of the B — > X s 7 given by the CLEO and 
Belle Collaborations [22] is a very powerful tool for constraining tan (3. We use the value [20] 

BR(6 -> si) = (3.37 ± 0.37 ± 0.34 ± 0.24+°;^ ± 0.38) • 10" 4 (21) 

where the errors are, respectively, the statistical error, the systematical error, the error from 
model dependence, the one due to the extrapolation from the data to take into account the 
full range of photon energies, and finally an estimate of the theory uncertainty. We allow the 
branching ratio to vary within the conservative range: 2 x 10 -4 < BR(6 — > 37) < 5 x 10" 4 . 

For the calculation of the MSSM Higgs boson and SUSY particle contributions, we will 
use the most up-to-date determination of the b — > S7 decay rate [36], where all known 
perturbative and non-perturbative effects are included 6 . This includes all the possibly large 
contributions which can occur at NLO, such as terms oc tan/5 [in particular in mj, and/or 
terms containing logarithms of M s /M w . For the input parameters we will use the values 
given in Refs. [36], except for the cut-off on the photon energy, E 1 > (1 — 5)m b /2 in the 
bremsstrahlung process b — > s^g, which we fix to 5 = 0.9. For the SUSY spectrum, we will 
work in the scenario discussed in section 2, i.e. a common squark mass M s = 1 TeV, the 
maximal mixing scenario At = v^Ms and A b = 1 TeV, fi = M 2 = 350 GeV. 



M A [GeV] 


90 


110 


130 


tan = 10 


5.2- 10" 4 


5.0- 10" 4 


4.8- 10~ 4 


tan p = 30 


3.7- 10" 4 


3.5- 10" 4 


3.3- 10" 4 


tan = 50 


2.5 • 10" 4 


2.3- 10~ 4 


2.2- 10~ 4 



Table 5: BR(b -> S7) in the MSSM for tan/3 = 10, 30, 50 and M A = 90, 110, 130 GeV. 

Values of BR(6 — > S7) in the MSSM are displayed for three choices of tan p and M A in 

Table 5, and one can see that they are in most cases within the allowed range 2 x 10~ 4 < 

BR(6 — > S7) < 5 x 10~ 4 . An exception is when tan/3 ~ 10 and the A boson is light, 

Ma ;$ 110 GeV, where the branching ratio slightly exceeds the maximal value; however, we 

have checked that for other values of the SUSY parameters a reasonable BR(b — > S7) can be 

accommodated also in this case. 

6 We thank the authors of Ref. [36], in particular P. Gambino, for providing us with their FORTRAN 
code for the calculation of BR(6 — > sj) at NLO. 
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4. Branching fractions and total decay widths 

The branching ratios and the total decay widths of the three neutral Higgs bosons [37], again 
evaluated with the program HDECAY with the set of input parameters discussed in section 
2.4, are shown in Figs. 4-5 and 6 as functions of their masses for Ma varying between 90 
and 130 GeV, for two values of tan (3, 10 and 30. 

In the case of the A boson, the situation is rather simple since in the entire mass range 
shown, the pattern is the same: BR(A — > bb) and BR(A — > t + t~) are approximately 
90% and 10%, respectively, while the gluonic branching fraction is slightly above 10~ 3 and 
BR(A — > ~ 5 x 10~ 4 . The branching ratios for the other decay modes, including 

the A — > 77 decay, are below the level of 10~ 5 for tan (3 > 10. This is due to the strong 
enhancement (suppression) of the A couplings to isospin — fermions [the Agg coupling 
is mediated here mainly by the 5-loop contribution]. Some of the final states which are 
important for larger Ma values [such as the decay A — > hZ] are strongly suppressed by 
phase space in addition to the coupling suppression for large tan (3 values. 

The pattern for the CP-even Higgs bosons is similar to that of the A boson if their masses 
are close to Ma and tan/? ^> 1. The exception is when h and H have masses very close 
to Mq. In this limit, they are SM-like and decay into charm quarks and gluons with rates 
similar to the one for t + t~ [~ a few %] and in the high mass range, M#-o ~ 125 GeV, into 
WW (~ 10%) and ZZ (a few %) pairs with one of the gauge bosons being virtual. However, 
this limit is not reached in practice, in particular for the lightest h boson with tan (3 < 50. 

Indeed, as can be seen in Fig. 4, all decays of the h boson other than bb and t + t~, are 
below the level of 10 -3 for tan f3 — 10 and 30. In particular, we have BR(/i — > 77) < 10 -4 
and BR(/i — > WW) < 10~ 2 even for M^ ~ 120 GeV, i.e. at least one order of magnitude 
smaller than in the case of the SM Higgs boson. This is a consequence of the fact that the 
h boson coupling to gauge bosons (and to up-type fermions) does not reach the SM limit 
while the couplings to 6-quarks are still enhanced even for Mh ~ Mq\ see Fig. 2 for the 
couplings. The gg branching fraction is in turn almost constant, BR(/i — > gg) ~ 10~ 3 , since 
the h — > gg decay width is dominantly mediated by the 6-quark loop, but is suppressed by 
a factor of more than 20 compared to the SM case. 

In the case of the heavier H boson, the trend depends strongly on the value of Ma- For 
M H > 125 GeV, the branching ratios into the 77, WW and gg final states are much smaller 
than for the SM Higgs boson: for tan (3 values larger than ~ 10, they are at least two orders 
of magnitude smaller. Near the critical mass M C) the branching fractions of the H boson are 
in general closer to those of the SM Higgs particle. This is particularly true for tan (3 values 
smaller than ~ 30, where the branching fractions into WW,gg and 77 (t + t~) final states 
are only slightly, less than a factor of two, smaller (larger) than for the SM Higgs particle. 
However, for larger values of tan/3 [not shown in Fig. 4a] a new feature occurs: the H 
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Figure 4: The branching ratios of the neutral Higgs bosons h, H and A as a function of their 
masses for two values of tanfi = 10 (left panel) and 30 (right panel). 
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boson coupling to down-type fermions which is proportional to tan (5 [in the approximation 
discussed in Section 2, g Hdd ~ —fi tan f3 + f 2 for Ma < M c ] becomes strongly suppressed 
and at some stage [close to tan/3 ~ /2//1 in our approximation] vanishes [note that this is 
not the case for up-type fermions and gauge bosons]. 

All these features can be seen more explicitly in Figs. 5, where the branching fractions 
of the h and H bosons into the t + t~, WW, gg and 77 final states are displayed as functions 
of tan/5 for three values of the A mass in the intense-coupling regime Ma = 90, 110 and 130 
GeV as well as in the decoupling limit with M4 = 1 TeV. In particular, one can see that 
the branching fractions BK(H — > t + t + ), and hence BR(iJ — > bb), drop at a given value of 
tan/3 [around tan/3 ~ 40 for M A = 90 GeV and tan/3 ~ 50 for Ma = 110 GeV], therefore 
enhancing the rates of the WW, gg and 77 decay modes. This is the only situation where 
the branching ratios for these decays are larger than in the case of the SM Higgs boson 7 . 

Finally, the total decay widths of the three neutral Higgs bosons are shown in Fig. 6. For 
h and H when they have masses close to Mc, the total decay widths are comparable to those 
of the SM Higgs boson, i.e. they are rather small. For tan (3 > 10, the A boson and one of 
the CP-even Higgs bosons, i.e. $a, decay most of the time into b quarks and r leptons and 
the total decay widths are to a good approximation given by: r tot ~ GpMA,<f> A /(4:y/2ir) x 
tan 2 f3 (to 2 + 3m 2 ) and are therefore rather large. For tan/3 = 50 and M A ~ 130 GeV, the 
total widths of the three neutral Higgs bosons are r tot (A) ~ 7.2 GeV, T tot {H) ~ 6.96 GeV 
(with M H = 127.6 GeV) and T tot (h) ~ 0.11 GeV (with M h = 123.9 GeV). 
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Figure 6: The total decay widths of the neutral Higgs bosons h, H and A as a function of 
their masses for two values of tanf3 = 10, 30. 



7 This "pathological" situation, where the couplings to down-type fermions are strongly suppressed, is 
known in the case of the the lightest h boson and is discussed in several places in the literature; see for 
instance Ref. [38]. 



24 



5. Production at Future Colliders 



In this section, we present the production cross sections of the neutral MSSM Higgs bosons 
in the intense-coupling regime at future colliders. We will consider the case of the hadron 
colliders, the Tevatron Run II and the LHC, a future e + e~ linear collider with a cm. energy 
of y/s = 500 GeV, in both the e + e~ and 77 modes [with the photons generated by back- 
scattering of laser light] and a /i + collider. In all cases, we will assume the values tan f3 — 10 
and 30 and a heavy superparticle spectrum with the soft SUSY-breaking scalar masses set 
to M s = 1 TeV and trilinear couplings A t = 2.6 TeV, A b — 1 TeV [we recall the reader 
that here, the higgsino and gaugino masses are set to /i — M 2 = 1 TeV]. Therefore only 
the standard processes, i.e. single or pair production and production in association with SM 
particles, will be considered. We will not discuss associated Higgs boson production with 
SUSY particles [39], Higgs boson production from the (cascade) decays of SUSY particles 
[40] or the effects of SUSY particle loops in the production processes [41]. 

5.1 Production at hadron colliders 

The main production mechanisms of neutral MSSM Higgs bosons at hadron colliders are the 
following processes [see Ref. [42] and for recent reviews Refs. [2, 43, 44] for instance]: 
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The single production of the pseudoscalar A boson in the weak boson fusion processes or in 
association with massive gauge bosons, as well as the pair production of two CP-even Higgs 
bosons, do not occur at leading order because of CP-invariance. Higgs boson production in 
association with other SM particles has too small cross sections. 

All the production cross sections are shown in Figs. 7-12 for the LHC with y/s = 14 TeV 
and for the upgraded Tevatron with y/s = 2 TeV. These cross sections have been obtained 
using the package CompHEP [45] [see Ref. [46] for details of the MSSM implementation in 
CompHEP by means of the LanHEP program], with the Higgs sector adapted from HDECAY 
for consistency with the discussion of the Higgs properties given previously. The parton 
densities are taken from the sets CTEQ5L and CTEQ5M1 [47] with the scale set at the 
Higgs boson mass. The choice of scale and parton density might lead to a 30% uncertainty 
in the production cross sections. 

The next -to-leading order QCD corrections are taken into account in the gg fusion 
processes where they are large, leading to an increase of the production cross sections by a 
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factor of up to two [48]. For the other processes, the QCD radiative corrections are relatively 
smaller [49]: for the associated production with gauge bosons and associated Higgs pair 
production, the corrections [which can be inferred from the Drell-Yan W/ Z production] are 
at the level of 10%, while in the case of the vector boson fusion processes, they are at the 
level of 30%. For the associated production with heavy quarks, the NLO corrections are only 
available in the case of the production of the SM Higgs boson with tt pairs [50] where they 
alter the cross section by ~ 20% if the scale is chosen properly. To include these corrections 
and to obtain the corresponding if -factors, we have linked CompHEP with the programs 
HIGLU, HPAIR, HQQ, V2HV and VV2H [51] which have all these corrections incorporated. 

Let us now briefly summarize the main features of these processes, relying for the exper- 
imental searches on the analyses performed by various collaborations [6-8]. 

a) Gluon gluon fusion 

This mechanism is mediated by top and bottom quark loops and is therefore sensitive to 
both the bb and tt-Higgs Yukawa couplings, allowing for the measurement of these impor- 
tant parameters [in conjunction with other production processes with the same Higgs decay 
modes] if the uncertainties from parton densities and higher order QCD corrections are prop- 
erly under control. In the large tan/3 regime, the 6-loop contributions are dominant since 
they are enhanced by tan 2 f3 factors. The top-quark loop gives a relatively important contri- 
bution only in the case of the SM-like CP-even Higgs boson i.e. h or H for Ma > Mc or 
Ma < Mc, respectively. The cross sections are not the same for the pseudoscalar A boson 
and the CP-even Higgs boson <$>a even for large tan f3 where the Higgs-66 Yukawa couplings 
are similar: this is due to the appearance of different loop induced Higgs-gg form factors for 
CP-even and CP-odd particles. The cross sections for the production of the A and h boson 
decrease quickly for increasing Ma [a factor ~ 5 for Ma varying from 90 to 130 GeV]; the 
cross section for H stays almost constant since Mu does not vary too much in this range. 

At the LHC, for tan f3 ~ 10, the cross sections for the A and h bosons vary from (9(150) 
pb for a low mass A boson, Ma ~ 90 GeV, to (9(30) pb for Ma ~ 130 GeV; the cross section 
for H production stays at the level of 30 to 50 pb in this mass range. This leads to a huge 
number of events, (9(10 7 ) per year for an integrated luminosity of / C ~ 100 fb _1 as expected 
in the high-luminosity option. For tan (3 ~ 30, the cross sections are one order of magnitude 
larger for the A and $a bosons [while it stays roughly constant for the boson] as a result 
of the increase of the Higgs-66 Yukawa coupling. At the Tevatron, the trend is practically 
the same as at the LHC [note that here, the initiated qq process plays a more important role 
than at high energies due to the reduced gluon luminosity] except that the cross section is 
~ 30 times smaller. Nevertheless, with the integrated luminosity of / C ~ 30 fb -1 expected 
at the end of the Run II option, this would lead to at least (9(10 4 ) Higgs bosons in each of 
the gg — > A, processes. 
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Figure 7: Production cross sections at the LHC (upper panel) and at the Tevatron (lower 
panel) for the gluon-gluon fusion mechanism as functions of Ma for tan/3 = 10, 30. 

Let us now discuss the signatures which would allow for the detection of these Higgs 
particles. In all cases, the bb final states cannot be used because of the huge QCD 2-jet 
background and the lack of a clean leptonic trigger. At the LHC, the final states which can 
be probed for the SM-like Higgs boson § H would be in principle the § H — > 77 final state 
for which the branching ratio can be of order 10~ 3 . However, this occurs only if is pure 
SM-like which, in the scenario discussed in Fig. 5, occurs only for the H boson for a rather 
light A boson and moderate to large tan (3 values. This is exemplified in Table 6 where are 
shown the cross sections for h, H and A boson production in the gg fusion mechanism times 
the Higgs decay branching ratio into 77 final states, relative to the SM case [with the SM 
Higgs boson having the same mass as the h, H or A boson]. In the case of H, only for large 
tan (3 and moderate Ma values, the cross section times branching ratio exceeds the SM value 
[with a factor up to 3]. For the pseudoscalar A boson, the very small 77 branching ratio is 
not compensated by the large cross section. Therefore, the ratio is much smaller than unity. 
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In the case of the lighter h boson, we are not yet in the decoupling limit even for M A ~ 130 
GeV and tan (3 ~ 50, and the cross section times branching ratio a(gg — > h) x BR(/i — > 77) 
is at least two orders of magnitude smaller than in the SM. 

Thus, there are situations where cr(<I>) x BR($ — > 77) is rather small for the three Higgs 
bosons h, H and A compared to the SM case, making this channel more difficult to use. In 
addition, for the A and $^ particles, the total decay widths are much larger than in the SM, 
leading to broader signals and therefore a much larger 77 background. The — > WW and 
ZZ final states would be also very difficult to use since in our case M<$, H ~ 130 GeV and the 
branching ratios are too small even for SM $ H VV couplings. 

At the Tevatron, although very difficult, only the — > t + t~ channel might be used as 
preliminary studies seem to indicate [52]. The branching ratios for the h, H and A boson 
decays into t + t~ final states are always of the order of 10% except in the "pathological" 
situations discussed in section 4, where the Ht + t~ coupling is suppressed. For the A and 
$a bosons, the channels A, <3>^ — * t + t~ can be used at the LHC in the high tan /3 regime, 
tan/3 > 10, where the cross sections are large enough. At the Tevatron, because of the 
smaller cross sections, this channel is expected to be more difficult, although the analyses 
mentioned previously give some hope. 



tan/3 


M A 


ct( j 4)xBR(77)|mssm 


M h 


ctO)xBR(77)|mssm 


M H 


o-(/T)xBR(77)|mssm 


f T(ifO)xBR(77)| SM 


-(i?0)xBR(77)| SM 


a(H°)xBR(77)| SM 




90 


1.5 • 10~ 3 


84.9 


7.5 • 10-3 


125.1 


0.6 


10 


110 


3.7- 10~ 4 


103.9 


5.6 • 10-3 


126.8 


0.15 




130 


1.4 • 10~ 4 


117.1 


1.2 • 10-2 


134.0 


1.4-10-2 




90 


1.8 • 10~ 2 


85.9 


2.2 • 10"2 


124.0 


1.9 


30 


110 


5.0 • 10~ 3 


106.6 


6.7-10-3 


124.2 


0.36 




130 


2.1 • IO- 3 


123.2 


1.2 • 10-2 


128.0 


1.8 • 10-3 




90 


5.2 • 10~ 2 


85.7 


5.6 ■ 10-2 


124.3 


2.0 


50 


110 


1.4-10- 2 


106.5 


1.5 • 10-2 


124.4 


3.2 




130 


6.2 • 10-3 


124.3 


2.3 • 10-2 


127.1 


4.3 • 10-3 



Table 6: The ratios a(gg — > $) x BR($ -> 77) in the MSSM with $ = h, H and A relative 
to the SM case for a SM Higgs boson with the same mass, for three values of tan (3 = 10, 30 
and 50 and three values of the pseudoscalar boson mass M4 = 90, 110 and 130 GeV. 
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b) Associated production with bottom quarks 

In these processes, the production cross sections for the A and <3>a bosons are also strongly 
enhanced by tan 2 (3 factors. The rates are similar to those in the gg fusion processes, (9(100) 
pb at the LHC and 0(1) pb at the Tevatron for M A ~ 90-100 GeV and tan (3= 10, and are 
one order of magnitude larger for tan/5 = 30. They decrease more slowly with increasing 
Higgs boson mass than in the case of gg fusion and the cross sections are closer in magnitude 
for the A and $,4 bosons [in fact, they must be equal for equal Yukawa couplings when the 
final state 6-quark mass is neglected, i.e. in the "chiral" limit]. In turn the cross section for 
associated production of the SM-like <3># boson is much smaller due to a tiny (not enough 
enhanced) Yukawa coupling. 
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Figure 8: Production cross sections at the LHC (upper panel) and at the Tevatron (lower 
panel) for Higgs bosons in association with b quarks as functions of Ma for tan/3 = 10, 30. 

Because of the small number of events, the detection of the boson in this process is 
practically hopeless. In turn, the detection of the A and $a bosons is more promising than 
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in the gg fusion process since we do not have to rely on decays into photons or massive gauge 
bosons which here, are absent or strongly suppressed. Indeed, the additional two 6-quarks 
in the final state, which can be rather efficiently tagged using micro-vertex detectors, will 
reduce dramatically the QCD backgrounds [especially since the production cross sections can 
be very large] to the level where the final states A, &a t + t~ at the LHC and A, §a bb 
at the Tevatron could be easily detectable for tan/3 > 10. 

c) Associated production with top quarks 
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Figure 9: Production cross sections at the LHC (upper panel) and at the Tevatron (lower 
panel) for Higgs bosons in association with top quarks as functions of Ma for tan f3 = 10, 30. 

Here, the production cross sections are suppressed by the smaller phase space [in particu- 
lar at the Tevatron] compared to the previous case and by the fact that the tt-Higgs Yukawa 
coupling is not enhanced compared to the SM case. To the contrary, the cross sections are 
strongly suppressed by tan 2 (3 factors for the production of the A and $a bosons. A reason- 
ably large cross section is obtained for the production of the boson, reaching the level of 
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~ 0.3 pb at the LHC and 2.5 fb at the Tevatron. In a large part of the parameter space we 
are concerned with, i.e. for 90 < Ma ^ 130 GeV and tan (3 > 10 [even around the turning 
point Ma ~ 125 GeV], the sum of the cross sections for h and H production is comparable 
to the one in the SM, since the sum of the squared couplings g\ t i + g 2 Ht i =1/ sin 2 (3 ~ 1 for 
tan p > 10 and the Higgs masses are comparable. 

The detection of the <& H boson in this process is in principle possible either through its 
77 or bb decay modes. In the 77 channel, the detection is however more difficult in general 
than in the SM, due to the reduced 77 branching ratio in most of the parameter space as 
discussed previously. The § H — > bb final state gives a more promising detection signal since 
the branching ratio is in general (slightly) larger than in the SM. 

d) WW/ZZ fusion processes 
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Figure 10: Production cross sections at the LHC (upper panel) and at the Tevatron (lower 
panel) for the vector boson fusion mechanisms as functions of Ma for tan/3 = 10, 30. 
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Here the cross sections are two to three orders of magnitude smaller than in the gg — > $ 
and bb& processes. They only reach the level of ~ 5 pb at the LHC and ~ 0.1 pb at the 
Tevatron for a SM-like Higgs boson which has maximal coupling to gauge bosons. For the 
$a boson, the rates are too small in particular for tan (3 > 10. 

At the Tevatron, the cross sections are very small and this process is too difficult to be 
used. At the LHC, the decay into t + t~ final states would allow for the detection of the 
particles, similarly as in the SM model, by taking advantage of the energetic quark jets in the 
forward and backward directions which allow for additional cuts to suppress the backgrounds. 
However, there are situations where both products a(VV h,H) x BR(h, H — > t + t~) are 
suppressed compared to the SM CclSC clS shown in Table 7 for several values of tan/3 and 
Ma- This occurs for large tan/5 > 30 and small Ma ~ 90 GeV, where the ratio <x($) x 
BR(tt)|mssm/o"(-^°) x BR(rr)| S M with $ = h(H) is strongly (slightly) suppressed. For 
tan/3 ~ 10 and Ma ~ 130 GeV, both ratios are slightly smaller than unity, 15% for $ = h 
and 30% for $ = H, but in this case the sum of the ratios is larger than one. In addition, in 
the pathological regions where the H boson couplings to r leptons is suppressed [for instance 
tan/3 ~ 50 and M A ~ 110 GeV], one can have a very tiny cross section times branching 
ratio, since in this case, the lighter h is $a like and does not couple to vector bosons. 
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Table 7: The ratios <j(VV -> $) x BR($ -> rr) in the MSSM with $ = h, H relative to the 
SM case for a SM Higgs boson H° with the same mass, for three values of tan /3 = 10,30 
and 50 and three values of the pseudoscalar boson mass Ma = 90, 110 and 130 GeV. 
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e) Associated production with W/Z bosons 

At the LHC, these associated production processes would allow for the measurement of 
the Higgs boson couplings to the massive gauge bosons [again in conjunction with other 
production and decay processes]. The cross sections for W$h associated production are 
at the level of ~ 1.5 pb at the LHC and ~ 150 fb at the Tevatron. They are of course 
slightly smaller than for a SM Higgs boson with the same mass. The cross sections for Z$ H 
associated production are smaller by approximately a factor of 2. Because of the fact that 
the gauge bosons are tagged through their leptonic decays, H^-boson final states are much 
more interesting since their leptonic branching fractions are larger. 




Figure 11: Cross sections at the LHC (upper panel) and at the Tevatron (lower panel) for 
the associated Higgs production with gauge bosons as functions of Ma for tan (3 = 10,30. 



At the Tevatron, these processes are the most promising ones to detect a SM-like Higgs 
boson at the Tevatron [44]. Since the branching fraction for the decay mode — > bb is 
always larger than in the SM, except in the pathological cases discussed above, the situation 
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is rather favorable for the process pp — > W§ H if the cross section is only slightly suppressed 
compared to the SM case. In fact, the ratio of cross sections times branching ratios into bb, 
a(QV) x BR(fofo), compared to the SM case, is the same as the one for the vector boson fusion 
with the Higgs bosons decaying into t + t~ final states [since the cross sections are governed by 
the same &VV couplings and the ratio BR($ — > r + r)/BR($ — > bb) ~ m^/3mf is constant]. 
Therefore, the situation is similar to that of the previous case (d) and the possibility that 
the cross sections times branching ratios are small for both h and H can occur. 

The situation is more complicated at the LHC if the Higgs bosons have to be detected 
through their 77 decay modes, since the 77 branching ratios are in general smaller than in 
the SM and the total decay widths of the states are much larger, as discussed previously 
[especially since already in the SM, the significance for this process is rather low]. 

f) Associated pair production 




90 100 110 120 130 90 100 110 120 130 



Figure 12: Cross sections at the LHC (upper panel) and at the Tevatron (lower panel) for 
the associated Higgs boson pair production processes as functions of Ma for tan (3 = 10, 30. 
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The associated pair production of CP-even and CP-odd Higgs bosons through the s- 
channel exchange of a Z boson has the lowest production cross section of all processes. The 
event rates are reasonable only for the qq — > Ah process in the lower mass range of the A 
boson, where the cross section is above the level of 0.1 pb at the LHC and 10 fb at the 
Tevatron. In this range, both the A and &a bosons will decay into bb pairs so that the most 
advantageous final state [in terms of production rate] to consider is 4 6-jets, the rate for the 
bbr + T~ final state being an order of magnitude lower. 

To our knowledge, the detection of MSSM Higgs bosons at the Tevatron and the LHC in 
these processes has not been subject to experimental scrutinity. This is probably due to the 
difficulty of extracting such a small 4 6— jet signal from the background. However, it might be 
interesting to look at such final states, not only for this particular process but also in Higgs 
boson pair production in higher order processes where the trilinear Higgs boson couplings 
are involved as will be discussed below. 

g) Pair production of Higgs bosons in higher order processes 

Pair production of Higgs bosons in the gg fusion occurs through triangle loop diagrams 
mediated by top and bottom quark exchange where a Higgs boson is produced (off-shell in 
our case) and splits into two Higgs particles, and box diagrams where the two Higgs particles 
are radiated from the internal top and bottom quark lines. These processes are of higher 
order in the electroweak coupling and therefore have small cross sections [53, 54] compared 
to the dominant processes discussed previously. However, they are interesting to consider 
since they involve the trilinear Higgs boson couplings [in the vertex diagrams where a virtual 
Higgs bosons splits into two Higgs bosons], the measurement of which is important in order 
to reconstruct the MSSM Higgs potential. We will thus discuss shortly these processes at 
the LHC and at the Tevatron. Additional processes such as the double Higgs production 
in association with W and Z bosons, the double Higgs production in the WW/ZZ fusion 
processes and the triple Higgs production in qq collisions, can also give access to the trilinear 
Higgs boson couplings; however, the cross sections are in general smaller than in the gg 
fusion mechanism and we will not discuss them here [54]. 

The production cross sections for the double production of Higgs bosons of any kind, 

(g) gg -> hh, Hh, HH, HA, hA, AA 

are shown in Figs. 13 for the LHC and Tevatron energies [for the latter case, only for 
tan/? = 30 and for those which are large enough]. At the LHC, the cross sections are in 
general well below the level of 0.1 pb. They can reach the picobarn level only for very large 
tan/5 values, tan/5 > 30, small A boson mass and only for the channels where only A or h 
are involved, gg — > hh and AA. In this regime, the Higgs bosons will decay predominantly 
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into bb pairs so that the final state will also be 4 6-quarks, as is the case for associated Higgs 
boson pair production which, as discussed previously, has similar cross sections in the same 
M A range and is a difficult process to detect experimentally. In fact, when the cross sections 
are sizable enough, the main contribution is due to the diagrams where both Higgs bosons 
are emitted from the internal 6-quark lines when the Higgs-66 Yukawa coupling is strongly 
enhanced. The triangle diagram which is sensitive to the trilinear Higgs boson coupling gives 
only a small contribution, and these couplings can thus not be probed. At the Tevatron, the 
cross sections are two orders of magnitude smaller than at the LHC because of the reduced 
phase space and the much lower gg luminosity. They are maximal for the gg — * hh and AA 
processes where they reach the 10 fb level in the same regime as at the LHC, i.e. for small 
Ma and large tan f5 where the contribution from Higgs emission from the internal 6-quarks 
is enhanced and the trilinear couplings play a minor role. 
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Figure 13: Production cross sections at the LHC (full lines) and the Tevatron (dashed lines) 
for the double Higgs boson production processes gg — > $i$2 as functions of Ma for tan (3 = 10 
(left panel) and tan/5 = 30 (right panel). 



5.2 Production at electron-positron colliders 

At e+e~ linear colliders operating in the 500 GeV energy range, the main production mech- 
anisms for MSSM neutral Higgs particles in the e+e- option are [2, 55, 56] 



(a) bremsstrahlung process 

associated production process 
(6) WW fusion process 

ZZ fusion process 
(c) radiation off top quarks 

radiation off bottom quarks 



e+e" -> (Z*) -> Z + h, H 

e+e- -> (Z*) -> A + h,H 

e + e- -> v e v e W*W* -> v e v e + h, H 

e+e- -> e + e'Z*Z* -> e+e" +h,H 

e+e- -> (7*,Z*) -+it_+h,H,A 

e+e~ -> -+bb + h,H,A 
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Again, the associated production of the pseudoscalar A and Z bosons, the single produc- 
tion of the A boson in the fusion processes as well as the pair production of two CP-even or 
two CP-odd Higgs particles can only occur at higher orders [57] because of CP invariance. 
The production cross sections are shown in Figs. 14-16 as functions of the pseudoscalar A 
boson mass for a center of mass energy ^/s = 500 GeV and for two values tan (3 — 10 and 
30. Let us briefly discuss the cross sections of these various processes. 

a) Higgs-strahlung and Higgs pair production 

These are the most interesting processes in this context. The cross sections for the 
bremsstrahlung and the pair production processes as well as the cross sections for the pro- 
duction of the h and H bosons are mutually complementary, coming either with a coefficient 
sin 2 (/5 — a) or cos 2 (f3 — a). The cross section for hZ production is large for large values of 
M h) i.e. close to the decoupling limit. By contrast, the cross section for HZ production is 
large for a light A boson implying a light H boson with a large coupling to Z bosons. For 
the associated production, the situation is opposite: the cross section for Ah is large for a 
light h boson whereas AH production is preferred in the complementary region. 
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Figure 14: Total cross sections at a 500 GeV e + e collider in the Higgs strahlung and 
associated Higgs production processes as functions of Ma for tan (3 = 10, 30. 

At y/s = 500 GeV, the sums of the cross sections <r(e + e~ — > hZ + HZ) and a(e + e~ — > 
hA + HA) are, respectively, ~ 60 and ~ 40 fb in the entire M A range. This means that with 
the integrated luminosity, f C ~ 500 ftr 1 , expected at the TESLA machine for instance, 
approximately 30.000 and 20.000 events per year can be collected in these two channels, 
respectively. In the Higgs-strahlung process, the signals consist mostly of a Z boson and a 
bb pair, which is easy to separate from the main background, e + e~ — > ZZ [for Mh — Mz, 
efficient b detection is needed]. For the associated Higgs boson production, the signals consist 
mostly of four b quarks in the final state, requiring efficient 6-quark tagging. 
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b) Vector boson fusion processes 

The trend for the production cross sections of the h and H boson in vector boson fusion is 
the same as in the Higgs-strahlung processes since both are proportional to the square of the 
Higgs-F V couplings. At y/s ~ 500 GeV, the production rate for the WW fusion mechanism 
is slightly larger than the one for the Higgs-strahlung process, a(e + e~ — > Huu + hvv) ~ 80 
fb compared to <x(e + e~ — > hZ + HZ) ~ 60 fb. A large sample of events is therefore expected, 
40.000 events per year with the luminosity f C ~ 500 fb -1 , with a signature consisting mainly 
of a bb or a t + t~ pair and missing energy. The cross sections for the ZZ fusion processes are 
one order of magnitude smaller than the ones for WW fusion, a consequence of the smaller 
neutral current couplings compared to charged current couplings. The signal is however 
cleaner due the additional e + e~ pair in the final state. 

Note that for other cm. energies, the Higgs-strahlung cross section, which scales as 1/s, 
dominates at lower energies while the WW fusion mechanism which has a cross section rising 
like \og(s/Mjj) dominates at higher energies. 



a(e*e~ —> X + h/H) [fb] 
Vs = 500 GeV 
tan^=10 




a{e*e" —> X + h/H) [fb] 
Vs = 500 GeV 
tan^ = 30 



Figure 15: Production cross sections at a 500 GeV e + e collider in the vector boson fusion 
processes as functions of Ma for tan/3 = 10, 30. 

c) Higgs production with heavy quarks 

For Higgs boson production in association with top quarks, the production cross sections 
are strongly suppressed for the A and $^4 boson for tan f3 > 10 and are sizeable only for the 
boson which has almost SM-like couplings to the top quarks. Even in this case, they 
are however very small, barely exceeding the level of 0.2 fb. This is due to the fact that at 
\fs = 500 GeV, there is only a little amount of phase-space [~ 50 GeV] available for this 
process. At higher energies, e.g. \fs = 800 GeV, the cross sections can exceed the level of 
~ 1 fb which would allow for the measurement of the $#tt couplings since most of the cross 
section is coming from the Higgs boson radiation off the top quarks. 
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In the case of Higgs boson production in association with bottom quarks, we have taken 
into account only the gauge invariant contribution coming from Higgs boson radiation off the 
6-quark lines since a much larger contribution would come from the associated production 
process e + e~ — > Ah or AH, with one of the Higgs boson decaying into 66, or from the Higgs- 
strahlung process e + e~ — > Zh or ZH with Z — > 66. These resonant processes have been 
discussed earlier and can be separated from the Higgs radiation off 6-quarks by demanding 
that the invariant mass of a 66 pair does not coincide with a Z boson or another Higgs boson. 
Because of the strong enhancement of the 66 Yukawa coupling, the cross sections can exceed 
the level of a(e + e~ — > bbA + 66$^) > 0.7 fb for tan/3 > 30, allowing for the possibility to 
measure directly tan/5 with a reasonable accuracy in this process [58]. 




d) Double Higgs production in the strahlung process 

Finally, we also show the cross sections at = 500 GeV and tan/? = 10,30 for the 
double Higgs boson production in the bremsstrahlung processes [59] 

(d) Double Higgs - strahlung e + e~ -> Zhh, ZHh, ZHH, ZAA 

which are the largest [the other production processes with more than one Higgs boson in 
the final state, such as WW fusion or triple Higgs production via Z exchange, have too low 
rates at this energy]. In particular, a(e + e~ — > Z$h&h) can be of the order of 0.1-0.2 fb 
[i.e. ~ 50 — 100 events for J C ~ 500 fb -1 ] which would allow a reasonable measurement 
of the & 3 H trilinear coupling with some accuracy [60]. The measurement of the other Higgs 
boson self-couplings is more challenging. 
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Figure 17: Production cross sections at a 500 GeV e + e~ collider in the double Higgs produc- 
tion processes e + e~ — > Z$i$ 2 as functions of Ma for tan (5 — 10 (left panel) and tan (3 = 30 
(right panel). 

5.3 Resonant Higgs production at 77 colliders 

77 colliders, in which the high-energy photon beams are generated by Compton back- 
scattering of laser light [61], provide useful instruments to search for the neutral Higgs 
bosons and to test their properties [62, 63]. Center of mass energies of the order of 80% of 
the e + e~ collider energy and integrated luminosities f C ~ 100 fb, as well as a high degree 
of longitudinal photon polarization can be reached at these colliders [64]. 

Tuning the maximum of the 77 spectrum to the value of the Higgs boson mass, which is 
assumed to be already known from the e + e _ option with some accuracy, the Higgs particles 
can be formed as s-channel resonances, 

77 -> h,H,A 

decaying mostly into bb pairs. The main background, 77 — > bb, can be suppressed by choosing 
proper helicities for the initial electron, positron and laser photons which maximizes the 
signal cross section, and eliminating the gluon radiation by taking into account only the 
two-jet cross section. 

Following the analysis of Ref. [63], to which we refer for details, we show in Fig. 18, 
the cross sections for resonant two-jet signal 77 — > h, H, A — > bb production. A cut on the 
scattering angle of the b quark has been applied, | cos#| < 0.5, which increases the ratio of 
the signal cross section to the 77 — > bb two-jet background cross section. [Note that the 
resummed NLO QCD corrections to both the signal, background and the interference have 
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been properly taken into account in the 2-jet final state]. The Compton spectrum has been 
integrated around the resonant Higgs mass in bins of A = ±3 GeV in order to account for 
the limited bb mass resolution. Clear signals can be obtained for tan (3 = 10 and 30, in the 
90 GeV < Ma ^ 130 GeV mass range, except for the pseudoscalar A boson for tan/5 = 10 
since in this case, the A77 coupling, mediated only by b and t quark loops, is suppressed. 




90 100 110 120 130 90 100 110 120 130 

Figure 18: Total cross sections for the resonant Higgs production 77 — > h, H, A — > bb, in 
Compton backs cattered 77 fusion for tan f3 — 10 and 30. The maximum of the 77 sub-energy 
is tuned to the mass of the $ boson. 

5.4 Resonant Higgs production at colliders 

The ability of a future muon collider to investigate the Higgs sector of the SM and MSSM 
has been discussed in numerous papers; see for instance Refs. [65-67] . The main advantages 
of a muon collider, compared to an electron-positron machine, are due to the fact that the 
muon has a much larger mass than the electron, which means that: 

(?) The couplings of Higgs bosons to pairs are much larger than the couplings 

to e + e~ pairs, yielding significantly larger rates for s-channel Higgs boson production at a 
muon collider [the production rate is negligible in e + e~ collisions]. 

(ii) A very narrow beam energy spread can be achieved, which leads to only a rather 
small loss of production rates at muon colliders. 

However, there is a dependence between the small beam energy spread, R = dE/E, 
and the delivered luminosity. Present estimates for yearly integrated luminosities are / C = 
0.1, 0.22, 1 fb _1 for beam energy resolutions of R = 0.003%, 0.01%, 0.1%, respectively [66]. 

Due to these features, a muon collider can be considered as a potential factory for the 
MSSM Higgs bosons, which can be produced as s-channel resonances: 

— > h, H, A 
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This would allow for the determination of the masses and the total decay widths of the 
MSSM neutral Higgs bosons with a high degree of accuracy. 

Following Ref. [67], we use a Gaussian spread for the center of mass energy, for 
which the central value is set at the Higgs boson mass, and calculate the effective s-channel 
production cross section for a Higgs boson <£> decaying into a final state X, using the formula: 

4tt BR($->/i+/QBR($->X) 



f(^/r tot ($ 



where y/s = M$, r to t( < ^) is the Higgs boson total decay width and is given by 



^ = aoo2GeV (aoS%)(lo^v) (23) 

In the discussed region of the MSSM parameter space for the intense-coupling regime, the 
three neutral Higgs bosons h, H and A are broad enough. Therefore, one can use a resolution 
R = 0.01% without too much loss of production rates. In such a case, the cross-sections are 
functions of the Higgs branching fractions and Higgs masses and practically do not depend 
on the resolution. The obtained cross sections as functions of the pseudoscalar Higgs boson 
mass Ma for tan/3 = 10 and 30 are shown in Fig. 19 [using eq. (22) with R = 0.01%]. As 
can be seen, they are large enough to allow for the production of a significant number of 
Higgs bosons, so that detailed studies of the profile of these particles can be performed. 




Figure 19: Cross sections for resonant MSSM Higgs boson (h,H and A) production at a 
H + H~ collider as functions of Ma for tan/5 = 10 (left panel) and 30 (right panel). 
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6. Conclusions 



In this paper, we have performed a comprehensive analysis of the MSSM Higgs sector in 
the intense-coupling regime, that is when the neutral Higgs particles are rather light, 90 
GeV < Mj, < 130 GeV, and the value of the ratio of the vevs of the two MSSM Higgs 
doublet fields, tan/3, is large. In this scenario, one of the CP-even Higgs particles, has 
SM-like Higgs properties, i.e. couples almost maximally to massive gauge bosons and top 
quark pairs, while the other CP-even Higgs particle, <!% has pseudoscalar-like properties, 
i.e. couples strongly to bottom quarks and r leptons and has almost no couplings to massive 
gauge boson pairs; however the ZQ^A coupling is almost maximal. 

We have first discussed the parameterization of the MSSM Higgs sector and studied the 
radiative corrections due to the top/stop and bottom/sbottom sectors to the Higgs boson 
masses, to their couplings to gauge bosons and fermions [as well as to their self-couplings]. 
We have derived rather simple expressions for these radiative corrections, which approxi- 
mate rather accurately the full corrections, and which allow to have a simple qualitative 
understanding of the behaviour of these various parameters in this regime. 

We have then analyzed the various experimental constraints on this intense-coupling 
regime, with large values of tan f3 and relatively small Higgs boson masses. We have shown 
that this scenario is still allowed by indirect constraints from high-precision measurements 
such as the electroweak parameters, i.e. the W boson mass, the effective electroweak mixing 
angle sin 2 6w and the Z boson partial width and forward-backward asymmetry in bb final 
states, the radiative and loop induced decay b — > 37 and the anomalous magnetic moment 
of the muon, (g — 2) M . We have also discussed the direct constraints from MSSM Higgs 
boson searches at LEP2 and the Tevatron. From Tevatron searches for the charged Higgs 
particle in top quark decays and the neutral Higgs bosons in the associated production with 
bb pairs, tan (5 values slightly above 50 are still allowed. Imposing the LEP2 constraints, we 
have delineated the allowed regions of the parameter space in tan/5 and M A (and M h ) in 
various scenarii for the mixing in the stop sector. We have also studied the implications of 
the possible ~ 2a evidence for a SM-like Higgs boson with a mass around 115 GeV, and 
have shown that in the intense-coupling regime, there are regions of the parameter space 
where this Higgs particle can be either the lighter h boson [for small to moderate mixing] or 
the heavier CP-even Higgs particle H [in the no-mixing scenario]. 

For the decay modes, the pattern for the CP-even Higgs particles is in general similar to 
the one for the pseudoscalar Higgs particle, i.e. they decay dominantly into bb [90% of the 
time] and t + t~ [10% of the time] pairs. There are three exceptions to this situation: 

(i) When the heavier H boson mass is close to its minimal value. In this case, it decays 
as the SM Higgs particle, i.e. the decays into cc, gg and WW* reach the few percent level. 
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(ii) When the lighter h boson mass is close to its maximal allowed value, the decays to 
other final states than bb and r + r~ are also sizeable. However, since the decoupling limit is 
not yet reached, the branching ratios are in general smaller than for a SM Higgs boson. 

(in) There are regions of the parameter space where the H boson couplings to down-type 
fermions are strongly suppressed, thus enhancing the decays into other final states. 

Thus there are regions of the parameter space where the decays of both the h and H 
bosons into the important channels 77 and WW* [and even r + r~ in some pathological 
situations] are suppressed compared to the SM case. The total decay widths of the A boson 
and one of the CP-even Higgs boson, $^4, are rather sizeable for large tan/5 values [a few 
GeV], while the total decay width of the other CP-even Higgs particle, is rather small 
as for the SM Higgs boson [of the order of 100 MeV]. 

A significant part of the paper was devoted to the analysis of the production of the three 
neutral MSSM Higgs bosons at the up-coming colliders: the Tevatron Run II, the LHC, a 
500 GeV electron-positron linear collider in both its e + e~ and 77 modes, as well as a future 
collider. All the neutral Higgs bosons [as well as the charged Higgs particle] will be 
kinematically accessible even at the upgraded Tevatron and at an e + e~ collider with a cm. 
energy ^Js > 300 GeV. The Higgs particles can be produced in various channels, allowing 
for a detailed study of the MSSM Higgs sector and for the determination of the various 
parameters such as the masses and the couplings. 

At e + e~ colliders, thanks to the clean environment, all Higgs particles can be easily 
detected and the Higgs boson masses and couplings can be measured with a relatively good 
precision. In particular, the couplings to photons, which could be sensitive to new particles, 
can be accurately measured in the 77 option of the collider. All production channels will 
be effective: the Higgs-strahlung and vector fusion processes for the SM-like & H boson 
and the associated Higgs pair production for the pseudoscalar A and pseudoscalar-like $a 
particles. In addition, both the associated production with top quarks [for the boson] 
and associated production with bottom quarks [for the A and $a particles] are possibly 
accessible, allowing for the measurement of the Higgs boson Yukawa couplings to top and 
bottom quarks and the determination of the important parameter tan/5. The total decay 
widths of two Higgs particles would be large enough to be resolved experimentally. The total 
decay widths of all three Higgs bosons can be precisely measured at a collider. 

At hadron colliders, the search for the MSSM Higgs particles in the intense-coupling 
regime can also be performed in different channels. At the Tevatron Run II with a very 
high luminosity, the SM-like <E># boson would be accessible in the associated pp — > 0>hV 
[V = W, Z] production process. In some cases, the production rates are larger than in 
the case of the SM Higgs boson; however, there are situations where the production rates 
are smaller for both the h and H bosons compared to the SM case [as is the case in the 
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pathological situation discussed for the H boson], making the detection more difficult. The 
$a and A particles can be produced in association with bb pairs with large rates for large 
enough tan (3 values and the states can be detected in most of the parameter space. 

At the LHC, all Higgs particles can be produced with rather large cross sections and many 
complementary production channels can also be effective. In particular, the production rates 
for the A and $^ bosons in the gg fusion mechanism and the associated production with bb 
pairs are strongly enhanced compared to the SM case, while the particles can be accessed 
in all the processes, gg, VV fusion and associated production with top quarks and vector 
bosons. However, the experimental detection of the particles can be slightly complicated for 
three reasons [this also holds in some cases for the upgraded Tevatron]: 

(i) Since the three neutral Higgs particles are relatively close in mass, it might be difficult 
to separate between them, in particular in the case of the A and $^4 bosons which mainly 
decay into bb and r + r~ final states which have a relatively poor resolution. 

(ii) The distinction between these two [A and $a] Higgs bosons can be made slightly 
more difficult by the fact that their total decay widths can be rather large, a few GeV, 
making the signal peaks rather broad. 

(Hi) The clean 77 (as well as WW*) final state signatures can be much less frequent 
than in the case of the SM Higgs particle, in particular when the lighter h boson plays the 
role of Again, the wider Higgs states would make the detection more difficult. 

In conclusion: the phenomenology of the MSSM neutral Higgs sector in the intense- 
coupling regime is extremely rich. All Higgs states are kinematically accessible at the next 
high-energy colliders, in various and complementary production processes. In some cases, 
the techniques for searching these particles will be different from the ones discussed in the 
case of the SM Higgs boson and even in the context of the MSSM Higgs bosons close to the 
decoupling limit. In this preliminary analysis, we have summarized the main features of these 
searches. However, detailed Monte-Carlo studies, including a more complete discussion of 
the various signals, a simulation of all the possible backgrounds, as well as a proper account 
of the experimental situation, will be needed to assess the potential of the colliders to discover 
these particles, and once discovered, to measure their properties. 
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